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Abstract 

We have found that some periods of interlayer exchange coupling (IEC) 
oscillations as a function of cap-layer (CL) thickness may be suppressed if 
the in-plane extremal spanning vectors of the cap- and ferromagnet-materials 
Fermi surfaces do not coincide. The suppression of the IEC oscillations vs. 
the CL thickness holds also if the magnetic slab thickness tends to infinity. On 
the one hand, we have shown by means of very simple arguments that apart 
from the well-known selection rules concerning the spacer- and cap-layers, 
another one related with the magnetic sublayers has to be fulfilled in order 
that the interlayer coupling oscillations vs. CL thickness could survive. On 
the other hand, the distribution of induced magnetic moments across the non- 
magnetic cap- and spacer-sublayers have been computed and shown to reveal 
the underlying periodicity of the materials they are made of (i.e. related to 
their bulk Fermi surfaces) independently of whether or not the selection rules 
are fulfilled. This means that the IEC oscillations are of global nature and 
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depend on all the sublayers the system consists of. 
75.70.i - Magnetic films and multilayers 

71.70Gm - Exchange interactions 

75.30Pd - Surface magnetism 
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I. INTRODUCTION 



Magnetic multilayers have been intensively studied for over a decade no wire. The reasons 
for it, apart from challenging cognitive aspects, are (already partially realized) practical 
applications of superlattices as magneto-resistive sensors, angular velocity meters, recording 
heads and magnetic memory elements. The phenomenon most of these applications is based 
on is the well known giant magneto-resistance (GMR) coming from a strong electron-spin 
dependence of resistivity in magnetic systems. To optimize devices of that sort, it is necessary 
to test the effect of all of the ingredients of the system in question (including kind of materials 
they are made from and thicknesses of particular sublayers), either directly on GMR or 
indirectly on the interlayer exchange coupling (IEC). Obviously, the effect of a spacer on 
IEC was established first&@, the next in turn was that due to magnetic sublayerdHH and 



finally the cap-layer (CL) effect has been studied quite recently^TEj. 

Before we present our original results let us briefly recall what are the most important 
facts concerning the CL's : i) the IEC oscillates as a function of CL thickness with a 



oscillations are strongly suppressed if stationary in-plane spanning vectors of the CL Fermi 



films shows a periodic distribution of the so-called quantum well states (QWS) with periods 
determined by extremal spanning vectors of the overlayer Fermi surface. We shall refer to 
the latter only indirectly, by exploiting the fact that the QWS lead to some spin-polarization 
of non-magnetic cap-layers. 

The aim of the present paper is to emphasize the relevance of magnetic sublayers to IEC 
oscillations as a function of cap-layer (CL) thickness. Besides, we shall comment on induced 
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magnetic moments in the non-magnetic sublayers, which may be viewed as a manifestation 
of the quantum well statesllUll. 



II. METHOD 

Our earlier paperlffi based on the single-band tight-binding model have proved that 
the model we use gives a reasonable qualitative description of basic physical mechanisms 
responsible for oscillatory phenomena in magnetic trilayers. 

Our Hamiltonian, described in detail in Ref. consists of the nearest-neighbor hop- 
ping and spin-dependent on-site potential terms. The systems under consideration now are 
trilayers capped with an overlayer, of the type n ovr O/rifF/n s S/nfF, where n ovr , rif and 
n s stand for the numbers of cap- (O), ferromagnet- (F) and spacer- (S) monolayers in the 
perpendicular z-direction. Hereafter the subscripts and superscripts ovr and s will always 
refer to the cap- and spacer-layers, whereas the spin- dependent parameters referring to fer- 
romagnetic sublayers will be indexed by o =f or |. For simplicity, we restrict ourselves to 
a simple cubic structure and regard the lattice constant and the hopping integral as the 
length- and energy-units, respectively. 

The interlayer exchange coupling has been calculated from the difference in thermo- 
dynamic potentials exactly as in Ref. ^l], moreover the magnetic moments (including the 
induced ones), m, have been expressed in terms of the eigen-f unctions u of the Hamiltonian 
as rrii = — n^, with n ia = \ u i a{E)\ 2 , where the summation runs over occupied states. 

E 

III. ASYMPTOTIC LIMITS 

In this section we present some analytic formulae which will be useful for interpreta- 



tion of rigorous numerical resluts of the next section. As has been shown in Ref. ^2|, the 
IEC can be Fourier-transformed with respect to n s and n a . That procedure can be quite 
straightforwardly generalized to include CL thickness as well. The resulting asymptotic 
(within the stationary phase approximation) expression consists of the terms of the form 



Apqrn\k\\i Ep)exp (2i(pk s z n s + (qk\ + rk^)nf + nk° vr n ovr )^j summed over all the in-plane wave 
vectors for which the exponential is stationary. The A coefficients are defined analogously 
as in Ref. |22|. Their exact numerical values are not important for qualitative considerations, 
one notes only that all the amplitudes of oscillations vanish asymptotically with the given 
sublayer thickness going to infinityil. There exist, however, some additional restrictions 
imposed by the asymptotic behavior of the IEC. In particular, a direct generalization of the 



results of Ref. g2j to the present case, with the cap-layer, gives: A 0qrn = (no coupling for 
n s — > oo). Another limit to be taken is rif — > oo, when, in view of the above mentioned 
asymptotic behavior, all the terms tend to zero except for A p q Q0 and A p00n . Since the os- 
cillations vs. spacer thickness survive in this limit in contrast to the ones vs. CL thickness 
which decay (see below), we conclude that A p000 ^ and A p00n = 0. 

Finally, taking into account the above mentioned restrictions and keeping for simplicity 
only the lowest order harmonics, we arrive at the following formula : 

J = ]T A W00 e 2ik ° n ° + A 1100 e 2i{Kns+krnf) + £ A 1010 e 2i ^ n ' +k ' n ^ 

a ai ct2 

_|_ A uoie 2i(kln s +kln f +k° vr n OV r) _|_ ^^ e 2i(k s -.n s +kin f +k° vr n ovr ) _|_ ^ 

where the a-s are the sets of in-plane wave vectors for which the relevant exponentials are 
stationary. For the case of the CL thickness dependence this allows us to formulate the 
following new selection rule, which in its general form (for n s and rif large and fixed and 
n ovr large and varying) reads : 

Vk° vr = 0, P Vk s z n s + (qVkl + rVk l z )n f = 0, (2) 

with nonvanishing p and either q or r (V is the two-dimensional gradient in the k x ,k y space). 
This means that out of all the stationary vectors of the cap material FS only those which 
simultaneously satisfy the above mentioned conditions for the in-plane gradients give rise 
to the oscillations with CL thickness. Eq. (0) is the main result of the present paper. 
This condition becomes even simpler in the particular case of the single-band simple cubic 
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model considered hereafter, when the second part of Eq. @ separates and all the individual 
in-plane gradients must vanish (c.f. Ref. |22|) . 

The origin of the new selection rule becomes clear if we qualitatively interpret Eq. ([]]) 
in terms of the quantum interference modefi The first term corresponds to the states 
reflected once at each of the spacer-ferromagnet interfaces, the second and third ones to 
the states penetrating one of the magnetic layers and reflected back at the cap-ferromagnet 
interface while the last two terms describe states reaching the outer boundary of the cap- 
layer ("vacuum"). It is quite clear therefore that the rij-dependent phase factor must be 
also taken into account while performing the stationary phase approximation. 

It is evident from formula (|l|) that the bias of oscillations with CL thickness depends not 
only on the spacer- and magnetic- layer thicknesses but on the on-site V ovr potential as well. 
The latter observation results from the fact that the A coefficients in the second and third 
terms of Eq. (|l|) depend on the value of the reflection coefficient at the cap-ferromagnet 
interface which in turn depends on the cap material electronic structure. 

The stationary spanning vectors, for a sublayer characterized by the potential V, can 
be determined in a very simple way, by minimizing with respect to k\\ the following Fermi 
surface equation for the s.c. lattice : 

k z (k\\, Ep) = arccos[(V — Ep)/2 — cos k x — cos k y ] . (3) 

Hence the in-plane extremal spanning vectors are: ku = (0, 0) for —6 < Ep — V < —2; 
(±tt, 0) (0, ±?r) for -2 < E F - V < 2; (±tt, ±tt) for 2 < E F - V < 6 and 

k z = arccosflV - E f ) /2 - a] , (4) 

with a = 2, and -2, for the corresponding kn, respectively. Thus, the period of oscillations 
vs. sublayer (with the potential V) thickness, is just A = ir/k z (or 7r/(7r — k z )). 
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IV. NUMERICAL RESULTS 



We shall now present our exact numerical results (see Ref. ^1] for details of the method) 
and show how they can be interpreted in terms of the analytical formulae from the preceding 
section. Fig. [l|. confirms the well-known fact that the IEC oscillations vs. CL thickness 
have got a period determined by kind of material the cap is made of and get suppressed if 
there is a mismatch in corresponding in-plane spanning vectors of the CL and the spacer. 
The suppression takes place in the cases c and d where k» — kn — (±7T, ±7r), opposite to 
fcover _ (o, ±7r),(±7r, 0). The dependence of the bias values on V ovr is also clearly visible. 
The new effect is presented in Fig. |], which shows that the suppression may be due to the 
misfit in the fen's corresponding to the overlayer and magnetic sublayers, respectively (curves 
c and d), whereas in case of the curves a and b the periodicity is quite pronounced owing 
to the matching of the above mentioned spanning vectors. It can be also readily seen from 
Fig. H] that the phases of oscillations as well as the bias-values depend on the potentials of 
the ferromagnetic layers (exchange splitting). It is noteworthy that Figs. [1] and|2| show that 
the selection rule works quite well, even when the relevant layer thicknesses are rather small: 
n s = 5 and nf — 10, respectively. This confirms our previous observation^ that relatively 
small systems in the z-direction may reveal the asymptotic behavior. A detailed inspection 
of curves c and d suggests that the selection rule is slightly more rigorously enforced in Fig. |2| 
(due to n/ = 10) than in Fig. |l| (due to n s = 5), but the effect is tiny indeed and hardly 
visible. Incidentally, all the periods of oscillations obtained by the numerical computations 
and visualized in Figs. |I]-|4] can be pretty well reproduced in terms of the asymptotic Eqs. @ 
and (|): e.g. for E F = 2.1 and V = -0.6, -0.3, 0.0, 0.3 and 0.6, we get A = 3.6, 4.9, 9.9, 6.9 
and 4.3 ML respectively. 

Another rather obvious but noteworthy effect consists in vanishing of the IEC oscillations 
vs. CL thickness when magnetic sublayer thickness gets bigger and bigger. This is shown 
in Fig. 0, and has not been discussed either, to our knowledge so far, although such a trend 



could be predicted on the basis of analytical formulae of Ref. [Tl. In fact this finding means 
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that in order to avoid undesirable effects of cap-layers (which may be of different thickness 
in an experiment) on the IEC oscillations one should work with thick magnetic sublayers. It 
is also noteworthy that the oscillation bias value depends on the magnetic layer thickness, 
as could be predicted from Eq. [I]. 

Finally, in connection with the quantum well states conceptEUll, we have studied the 
distribution of induced magnetic moments in the CL (and in the spacer). A typical result 
is presented in Fig. |j. The induced magnetic moments are measured in dimensionless units 
(fiB = 1) and are of the order of 0.1% with respect to the magnetic layer magnetization. 
As expected, the period of the induced-moment distribution within the CL is exactly that 
anticipated for the bulk CL material FS. The effect of the other sublayers is minor, except 
that the magnitude of the induced moments is also magnetic-slab dependent. This might 
seem, at a first glance, to be in conflict with the IEC behavior which shows no oscillations for 
parameters of Fig. [| (cf. Fig. |I]c). Yet the spin polarization in non-magnetic layers is related 
to just one system with the fixed sublayer thicknesses and the given alignment of magnetic 
sublayers, whereas the IEC results from the total energy (thermodynamic potential) balance 
between the two possible ferromagnetic layer alignments and has to do with the series of 
samples with changing CL thicknesses. This observation implies that the induced-magnetic 
moments in the non-magnetic cap-layer (as well as the QWS) give in general the whole set 
of periods, out of which only those survive, as far as the IEC is concerned, which fulfill 
the selection rules referring to the entire system. In other words the IEC oscillations are 
the global characteristic of the whole system, whereas the induced spin polarization in the 
cap-layer is strictly of local nature. 

The selection rules completed hereby by the extra condition related with the extremal 
spanning vectors of the magnetic sublayers, are quite general and apply to real systems, 
too. In particular they allow to explain why in case of the Cu/Co/Cu/Co multilayer the 
short period of oscillations with Cu cap-layer thickness is absent^ in spite of theoretical 
predictions^ and the photoemission results concerning QWS00E1. In fact the explanation 



is simple and quite analogous to that of Ref. |22] about IEC oscillations as a function of 



ferromagnetic layer thickness. Out of two in-plane extremal spanning vectors of the Cu 
Fermi surface only the "belly" one (at k» — ) coincides with the extrema of the majority 
and minority sheets of the Co Fermi surface, giving rise to the long period of oscillations. 
The "neck" spanning vector has no counterpart in the Co FS and this is why there are no 
short period oscillations. 

In conclusion, we have shown that in order for the interlayer exchange coupling os- 
cillations vs. cap-layer thickness to exist, it is necessary, that both the cap-layer- and 
magnetic-layer- Fermi surfaces share the same extremal in-plane spanning vectors. If this 
new "selection rule" is not fulfilled the period anticipated from the bulk cap-layer material 
will not occur in the exchange coupling, although it will still be present in the induced 
moment distribution across the cap layer. Another finding of this paper is that the IEC 
oscillations vs. CL - thickness vanish if magnetic sublayers thickness tends to infinity. 

This work has been carried out under the KBN grants No. 2P03B 165 10 (MZ) and 
2-P03B-099-11 (SK). We thank the Poznan Supercomputing and Networking Center for the 
computing time. 
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List of Figures 

1 The exchange coupling vs. cap layer thickness for n s = 5, rif = 3, Ep = 2.1. 
V s = = Q and V = —2.0. Stationary in-plane spanning vectors of the 



spacer and both the ferromagnetic FS's are k\\ = (±7T, ±7t). For the curves~a 
and b the stationary in-plane vector of the cap FS is the same in contrast with" 



the c and d cases, for which k\\ = (0, ±7r), (±7t, 0) what results in suppressing 
the oscillations. 



I IEC vs. CE thickness for = 5, w f = 10, A' F = 2.1, P = = O.d 



and = —2.0. Stationary in-plane spanning vectors for the spacer and the 
cap-layer are k \ \ = (0, ±7r), (±7r, 0). For the first two curves (a and 6) the 



minority spin FS stationary points coincide with those of the spacer and the 



overlayer. In case of the curves c and d both the majority and minority spin 
Fermi surfaces have the k\\ = (±7T, ±7r) spanning vector and consequently trie 
IEC oscillations are suppressed. 
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The effect of magnetic sublayer thickness on 


the IEC oscillations as a function 












of cap-layer thickness (parameters as in Fig 


. 1 except for V ovr = -0.6). 
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Induced magnetic moments (with /ig = 1) 


for n s = n om . = 40 and rif 


= 10 



for parallel (full line) and antiparallel (dashed line) configurations. The other 



parameters as in Fig, lc. Thick vertical lines mark the interfaces 
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FIG. 1. The exchange coupling vs. cap layer thickness for n s = 5, rif = 3, Ep = 2.1, 
V s = = and = —2.0. Stationary in-plane spanning vectors of the spacer and both 
the ferromagnetic FS's are fey = (±7r,±7r). For the curves a and b the stationary in-plane vector 
of the cap FS is the same in contrast with the c and d cases, for which fey = (0, ±7r), (±7r, 0) what 
results in suppressing the oscillations. 
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FIG. 2. IEC vs. CL thickness for n s = 5, n f = 10, E F = 2.1, V s = V ovr = 0.3 and V T = -2.0. 
Stationary in-plane spanning vectors for the spacer and the cap-layer are fen = (0,±7r), (±7r, 0). 
For the first two curves (a and b) the minority spin FS stationary points coincide with those of 
the spacer and the overlayer. In case of the curves c and d both the majority and minority spin 
Fermi surfaces have the k\\ = (±7r, ±7r) spanning vector and consequently the IEC oscillations are 
suppressed. 
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FIG. 3. The effect of magnetic sublayer thickness on the IEC oscillations as a function of 
cap-layer thickness (parameters as in Fig. 1 except for V ovr = —0.6). 
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FIG. 4. Induced magnetic moments (with //g = 1) for n s = n ovr = 40 and n/ = 10 for parallel 
(full line) and antiparallel (dashed line) configurations. The other parameters as in Fig. lc. Thick 
vertical lines mark the interfaces. 
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